Entanglement oscillations in non-Markovian quantum channels 
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We study the non-Markovian dynamics of a two-mode bosonic system interacting with two uncorrelated 
thermal bosonic reservoirs. We present the solution to the exact microscopic Master equation in terms of the 
quantum characteristic function and study in details the dynamics of entanglement for bipartite Gaussian states. 
In particular, we analyze the effects of short-time system-reservoir correlations on the separability thresholds 
and show that the relevant parameter is the reservoir spectral density. If the frequencies of the involved modes 
are within the reservoir spectral density entanglement persists for a longer time than in a Markovian channel. 
On the other hand, when the reservoir spectrum is out of resonance short-time correlations lead to a faster 
decoherence and to the appearance of entanglement oscillations. 
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The loss of coherence, or decoherence, within a quantum 
system is due to the interactions of the system with the sur- 
rounding environment. Decoherence processes cannot be ig- 
nored, for they are the main obstacle to the full realization 
of quantum information processing. The dynamics of open 
quantum systems, however, may be rather involved, mostly 
due to the complex structure of the environment interacting 
with the quantum system. Therefore, in order to describe the 
dynamics of the system of interest, some approximations are 
often made, leading to the derivation of a Master equation 
for the reduced density matrix. The most relevant approxi- 
mations are the weak coupling or Born approximation dyD, 
assuming that the coupling between the system and the reser- 
voir is small enough to justify a perturbative approach, and the 
Markov approximation |l|,|2|], which amounts to neglect short- 
time correlations between the system and the reservoir. There 
are few cases, however, where an exact analytic description 
of the dynamics is possible. Two relevant examples are the 
quantum Brownian motion (QBM) d |^, [H 01 and the case 
of a two-level atom interacting with a thermal reservoir with 
Lorentzian spectral density (ll |5| • 

Entanglement in continuous variable (CV) quantum chan- 
nels has attracted much interest in recent years, due to po- 
tential improvement in the channel capacity [6]. A realistic 
analysis of CV channels must take into account decoherence 
and dissipation phenomena and, in fact, there has been an 
increasing interest in the description of noisy CV quantum 
channels 0]. Most of the theoretical descriptions, however, 
rely on Born and/or Markov approximations and only very 
recently some phenomenological models of non-Markovian 
quantum channels have been proposed H,^. Non-Markovian 
effects are crucial, e.g., for high-speed quantum communica- 
tion where the characteristic time scales become comparable 
with the reservoir correlation time. Moreover, when the sys- 
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tem interacts with a structured reservoir, e.g., for quantum 
channels embedded in solid-state devices, memory effects are 
typically non negligible. In these cases the dynamics can be 
substantially different from the Markovian one. 



In this communication we focus attention on the dynam- 
ics of a two-mode bosonic quantum system propagating in 
a noisy bosonic quantum channel, i.e., interacting with two 
uncorrelated bosonic thermal reservoir. Such a system may 
be seen as made by two quantum Brownian oscillators and 
therefore it is possible to exactly describe the system dynam- 
ics. Using this analogy, we generalize the single mode QBM 
solution to the bimodal case and analyze in details the en- 
tanglement dynamics of bipartite Gaussian state. We present 
an exact approach based on a microscopic model and, start- 
ing from the exact solution, we discuss in details the genuine 
non-Markovian effects on the entanglement dynamics for the 
relevant class of bipartite Gaussian states. In particular we 
analyze the effects of short-time system-reservoir correlations 
on the separability thresholds. The exact results are compared 
with those obtained within the Markovian approximation, and 
few relevant situations in which the system-reservoir correla- 
tions give rise, for short times, to entanglement oscillations are 
addressed. We stress that our approach describes the effects of 
quantum noise for two-mode quantum systems starting from 
a microscopic description of the system and the reservoir. 



The exact Master equation describing a quantum harmonic 
oscillator interacting with a bosonic reservoir in thermal equi- 
librium has been derived for the first time in Ref. J3l and it 
is usually referred to as Hu-Paz-Zhang Master equation. The 
original derivation is based on the influence functional-path 
integral formalism though a simpler derivation based on the 
time-convolutionless projection operator technique (U may 
be also used [see, e.g., Ref. [10]]. Assuming that the bi- 
modal field interacts bilinearly with two identical uncorre- 
lated bosonic thermal reservoirs the Master equation for the 



2 



reduced density matrix g(t) of the field is given by 

=E{^°>e(*)]- A(t)[x*, [**,<?(<)]] 



n(<)[X fc ,[P fc ,^)]] + -r(t)[X fc 2 ,e(i)] 



i 7 (i)[Xfc,{P fc ,f>}]}, 



(1) 



where 



1 

3 



Pa. 



71 H 



a* 



are the dimensionless quadrature operators, [ofc,fflt] = 1 
(k = 1,2) being the mode operators of the two oscillators, 
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fiwoCofcafe + 1/2). This Master equation, being ex- 
act, describes also the non-Markovian system-reservoir cor- 
relations due to the finite correlation time of the reservoir. 
In contrast to other non-Markovian dynamical models fiUl . 
Eq. ([TJi is local in time, i.e. it does not contain memory inte- 
grals. All the non-Markovian character of the system is con- 
tained in the time dependent coefficients, A(t), O(t), r(t) and 
7(t), appearing in the Master equation. These coefficients de- 
pend only on the reservoir spectral density, i.e. on the micro- 
scopic effective coupling strength between the system oscilla- 
tor and the oscillators of the reservoir. The coefficient r(t) de- 
scribes a time dependent frequency shift, j(t) is the damping 
coefficient, A(t) and H(t) are the normal and the anomalous 
diffusion coefficients, respectively 01 [H]. It is worth under- 
lining that the Master equation given by Eq. (Q~|), is valid for 
general forms of the reservoir spectral density J(oj) and any 
temperature T, 

The solution of Eq. (fl} can be obtained generalizing to two- 
modes the method of solution of the Hu-Paz -Zhang Master 
equation. Several approaches have been proposed in the liter- 
ature to solve the Hu-Paz-Zhang Master equation fl2i[T3i[T4ll . 
The one developed in Refs. lfl3l [Till is based on the symmet- 
rically ordered quantum characteristic function (QCF) Xt(0 
at time t. Extending this approach to the case of a bimodal 
field, in the interaction picture with respect to J2k=i 2 we 
obtain the following solution 



X t(A) = cxp {-A T A t A} Xo ( e - r ^A), 



(2) 



where we indicate with xo(A) the characteristic function at 
t = 0, and where A = (xi,yi,X2,y2) T , (• • ■ ) T being the 
transposition operation. The bipartite QCF is defined as 



X(A) = Tr [ e «i a l-«i* Q i e « 2a 2-«2t 



Q 



(3) 



with £fc = -i=(xfe + iyk), and k — 1,2. From the QCF one 
evaluates the moments of the field 

{a\ k 4a™a%) = {-) m+n d^d^d^td^e^ 2 + ^ 

x x(6,6)ki=o,«2=o , 



and, in turn, the time evolution of the covariance matrix, thus 
quantifying the entanglement between the two modes. In 
Eq. (ffji the matrix A t is given by iflill 



-r(t) 



f ds e r(s) J R T (i, s) M(s) R(t, s), (4) 
Jo 



where the matrix R(t, s) contains rapidly oscillating terms 
and M(t) is given by 



M(t) 



I f 2A(t) -n(i) 

2 1 -n(t) o 



(5) 



It is worth recalling that the time dependent coefficients A(t), 
II(i) and j(t) depend only on the spectral density of the 
reservoir and can be expressed as power series in the system- 
reservoir coupling constant a. Finally, T(t) takes the form 



r(*) = 2 f ds 7 (s) 



(6) 



Once the spectral characteristics of the environment and the 
initial state of the field have been specified, then Eq. © fully 
characterizes the dynamics of the system without approxima- 
tions. 

We now consider a specific example of interest for CV 
quantum information and communication, namely the case of 
a bipartite Gaussian state of the form 



Xo(A) = cxp <{ -iA T er A - iA T X t 



where 



er = 
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C 
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(7) 
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with A = Diag(a, a), B = Diag(6, b), a,b > 0, and 
C = Diag(ci,C2) is the covariance matrix in the canonical 
form and X in = Tr[g(t) (Xi, Pi, X 2 , P 2 ) T ] the vector of 
mean values. For optical bimodal fields the solution of the 
Master equation can be simplified invoking the secular ap- 
proximation. It has been shown that, in this case, the time 
evolution depends only on the diffusion coefficient A(t) and 
on the dissipation coefficient 7 (i) ifTolfTsifTill . The expression 
for the QCF at time t becomes 



Xt(A) = cxp 



— -A ctjA — iA X t 



(9) 



where the evolved mean value and covariance matrix are given 
by 

X t = e- r ^/ 2 (R®R) T X ia , (10) 
<r t = e- rw (R®R) T (T a (R®R) + 5A r (*)l , (11) 
1 being the 4x4 identity matrix and R the rotation matrix 



R 



cos ojQt sin ojQt 
— sin ujot cos ujot 



(12) 
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Finally, the time dependent coefficient Ar(t) appearing in 
Eq. dTTb is given by 



A r (i) = e" r(t) / dse r(s) A(s) 



(13) 



The map expressed by Eqs. ( fTOb and ( fTTT i is Gaussian, i.e. 
it maintains the Gaussian character of the input state. As a 
consequence, separability after interactions may be checked 
by positivity of the partially transposed density matrix (PPT 
condition) fm[l7ll . which, in terms of the covariance matrix, 
reads as follows 



where f2 = u> © u) T ' with 




(14) 



(15) 



Explicitly, for the evolved state of Eq. (O, the PPT condition 
may be written as 



(a + 6)e" r W +A r (t)] 

-e- r «{[(a-&) a + 2(c? + c&] e- r W 
-2 ^/(a - 6) 2 + (ci - c 2 ) 2 + (c 2 - C 2)2 e -2r(t) | > 1 



(16) 



In the case of a twin-beam (TWB) state of radiation, i.e.. 



1 



c 2 = — sinh(2r) 



a = b = — cosh(2r) C\ 



r > 0, the last inequality reduces to [18] 

S(t) = e- 2r e- r W + A r (t) - 1 > 0. (17) 

Notice that, since in the case of a Gaussian state there is no 
bound entanglement, the sign of S(t) is strictly connected to 
the separability of the state: when S(t) > the state is sepa- 
rable and, for this reason, S(t) is called separability function. 
At t = 0, T(0) = A r (0) = 0, hence 5(0) = e~ 2r - 1, i.e. the 
initial state is an entangled state. As time passes, due to the 
interaction with the external environment, the entanglement 
between the modes deteriorates and the separability function 
S(t) vanishes. The time t s at which S(t s ) = 0, i.e. when the 
entanglement between the two modes is lost, is referred to as 
the separability threshold or separability time. Of course, t s is 
a function of both the system and the reservoir parameters. In 
the Markovian approximation, for example, t s has been shown 
to become infinite for a zero T reservoir II 1911 . 

In order to study in more detail the entanglement dynamics 
in a non-Markovian channel, we need to specify the spectral 
density of the reservoir. We focus on the Ohmic reservoir with 
Lorentz-Drude regularization, whose spectral density is ||2|] 



where uj c is the cutoff frequency. This is one of the most stud- 
ied models of bosonic environments since it leads to a friction 
force proportional to velocity, which is typical of dissipative 
systems in several physical contexts. In this model, the reser- 
voir correlation time is given by the inverse of the cutoff fre- 
quency, i.e. tji = l/ui c . For times t < tr the system dy- 
namics cannot be approximated by the Markovian dynamics, 
as we will show in the following. The expressions for the rel- 
evant time dependent coefficients, up to second order in the 
system-reservoir coupling constant, are given by ll4l. [l5ll 



A(t) = / dr k(t) cos(ujqt) 
Jo 

7(t) = / dr/i(r) sin(w r) 
Jo 



(18) 
(19) 



with 



k(t) = a 2 / duj J(uj)coth[huj/2k B T}cos(ujT) (20) 
Jo 

(21) 



/i(t) = a dojJ(uj) sin(wr) 



being the noise and the dissipation kernels, respectively |Q]]. 
We note that, for high reservoir temperatures T, A(i) ^> j(t). 
Inserting the expression of J(u>) into Eqs. ( fT8l and ( fT9l one 
gets analytic expressions for both of the coefficients for the 
high T case. For generic temperatures only j(t), which is in- 
dependent of temperature, has a simple analytic form, while 
the time integral in A(t) has to be calculated numerically 
[see, e.g., Ref. [15] for the analytic expression of the dif- 
fusion and dissipation coefficients]. We begin our analysis 
by noting that, for times t 3> tr, the coefficients A(£) and 
7(t) can be approximated by their Markovian stationary val- 
ues Am = A(t — > oo ) and jm — j(t — > oo) [see Appendix 
B of Ref. lfl5h l. Eqs. © and ( fl3l , hence, become 



r(t) = J M t 
A r (t) = [2N(lu ) 



1](1- 



-7m* 



(22) 
(23) 



where N(u>o) — ( e hu o/k B T _ ^-l j s ^ avera g e number of 
thermal photons. Inserting Eqs. d22l > and (f23b into Eq. ( TTTb . 
one can easily see that the separability time coincide with that 
predicted by the Markovian theory jl8lll9ll . i.e 



■ log 1 , 



Notice that, as mentioned above, if N(u>q) — > then t s — > oo. 

We now look at the deviation from the Markovian value of 
the separability time for times t < tr. In this time interval 
the expression for the separability condition ( TTTb reduces to 



S(t) 



-2r 



ds "f(s] 



dsA(s) - 1 > 0. 



7r uj;. + LU Z 



(24) 
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We consider the high T case. Using the analytic expressions 
of A(t) and j(t), we get 



S(r)=a — , , 



x 2 -l 



T ■ 



frlU c 1 + X 2 1 X 2 + 1 

a; 



1 — e r cos(t/ie) 



_ , e- T sin(r/x) +e"*-l, (25) 
x 2 + 1 J 

with t = w c i, and a; = uj c /ujq. In Fig.[TJwe report the separa- 
bility function S(t), and we compare it with the corresponding 
Markovian approximation 



Sm(t) 



21 + x 2 



huj r 



1 . (26) 
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FIG. 1: Separability function S(t) for twin-beam state in Markovian 
and non-Markovian channels for a high T reservoir. Plots refer to pa- 
rameters ksT/HuJc = 100, a 2 = 0.01, r = 0.1, and (Top): x = 10; 
(Bottom): x — 0.01. The top plot refers to a situation in which the 
frequencies of the involved modes overlap with the reservoir spectral 
density (x S> 1). Here , the more accurate non-Markovian esti- 
mation predicts that the entanglement persists for a longer time. This 
conclusion holds for any value of the reservoir temperature, provided 
the high T condition is satisfied, and as long as x 3> 1. The bottom 
plot refers to a situation with the reservoir spectrum out of resonance 
with respect to the mode frequencies, i.e. when x <C 1. Here, the 
environment causes a faster loss of entanglement compared to the 
Markovian prediction, as well as to the appearance of entanglement 
oscillations. 



The case considered in the top plot corresponds to a situa- 
tion in which the frequencies of the involved modes overlap 
with the reservoir spectral density (x 1). The Marko- 
vian prediction for the separability time, in this case, gives 
a smaller value than the one predicted by our non-Markovian 
theory. More precisely we have r s ~ 0.18 < r s ~ 0.65. 
Therefore, the more accurate non-Markovian estimation pre- 
dicts that the entanglement persists for a longer time. It is easy 
to show that this conclusion holds for any value of the reser- 
voir temperature, provided the high T condition is satisfied, 
and as long as x ^S> 1. 

The situation changes for the case in which the reservoir 
spectrum is "out of resonance" or "detuned" with respect to 
the frequency of the TWB, i.e. when x <C 1. As shown in 
Fig. 1 (bottom plot), in this case the environment causes a 
faster loss of entanglement compared to the Markovian pre- 
diction. A careful analysis of Eq. ( fT7b shows that this fea- 
ture can be traced back to the short time behavior of the dif- 
fusion coefficient A(i), describing environment induced de- 
coherence u3l. For x 3> 1, indeed, A(t) > Am, with 
Am = A(i — > oo) the Markovian value. Therefore, in 
this case, the system experiences a weaker decoherence than 
the one predicted by the Markovian approximation and, as a 
consequence, the entanglement between the modes persists 
longer. On the contrary, for x -C 1, A(t) quickly assumes 
values higher than Am- The system hence is subjected to a 
stronger initial decoherence which quickly destroys the quan- 
tum entanglement between the two modes. In Fig. 1 we see 
that, for x <C 1, the separability function displays entan- 
glement oscillations, a typical non-Markovian feature. Sim- 
ilar oscillations have been extensively studied for the case of 
the one-mode harmonic quantum Brownian model lf2Tj|l . It 
is worth noting, however, that here we look at quantum cor- 
relations between two modes whose dynamics cannot be di- 
rectly deduced from the dissipative dynamics of a single mode 
bosonic quantum channel. 

In conclusion, we have analyzed the dynamics of bipartite 
Gaussian states in a non-Markovian noisy channel. We found 
that when the mode frequencies are within the reservoir spec- 
tral density entanglement persists for a longer time than in a 
Markovian channel. On the other hand, when the reservoir 
spectrum is out of resonance than short-time correlations lead 
to a faster decoherence and to the appearance of entanglement 
oscillations. Our analysis indicates that non Markovian fea- 
tures of system-reservoir interaction may be relevant in quan- 
tum information processing and should be taken into account 
in the design of any quantum communication channel. 
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